Equilibration of Isolated Macroscopic Quantum Systems under Experimentally 

Realistic Conditions 
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In how far does an non-equilibrium initial ensemble evolve towards a stationary long time behavior 
for an isolated macroscopic quantum system? We demonstrate that deviations from a steady state 
indeed become unmeasurably small or exceedingly rare after initial transients have died out under 
the following conditions: The Hamilonian does not exhibit exceedingly large degeneracies of energy 
eigenvalues and energy gaps. A large number of energy levels is significantly populated by the initial 
ensemble. The system is observed by a measurement device with a reasonably bound working range 
compared to the resolution limit. The entire experiment, ending with the quantum mechanical 
measurement process, can only be repeated a "reasonable" number of times. It is argued that 
all these prerequisites for equilibration are fulfilled under many, if not all, experimentally realistic 
conditions. 

PACS numbers: 05.30.-d, 05.30.Ch, 03.65.-w 
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I. INTRODUCTION 

According to textbook Statistical Physics, all prop- 
erties of a macroscopic system at thermal equilibrium 
are perfectly described by the canonical ensemble when 
weakly coupled to a thermal bath and by the microcanon- 
ical ensemble when isolated from the rest of the world. 
However, there still seems to be no truly satisfying ex- 
planation of why this is so [l|, |2[ . 

Here we adopt the most common viewpoint that in- 
vestigations of this problem should be based on standard 
Quantum Mechanics and should start with the treatment 
of isolated systems. Indeed, it is widely believed that by 
considering a system in contact with a thermal bath as 
a single, isolated "supersystem" , the canonical formal- 
ism should be deducible from the microcanonical one. In 
turn, to directly describe an open system alone (without 
the bath to which it is coupled) by means of standard 
Quantum Mechanics does not seem possible. 

Excellent reviews of the enormous number of pertinent 
"classical" (i.e. older) works are provided e.g. by [H, 
while more recent developments are covered e.g. by 
With our present work we summarize and further develop 
a recently developed new approach frjl-fl2j to the funda- 
mental question: In how far does an non-equilibrium ini- 
tial ensemble evolve towards a stationary long time be- 
havior? Besides this problem of equilibration, a second 
fundamental issue is the problem of thermalization, i.e., 
the question whether, and to what extend, such a steady 
state is in agreement with the microcanonical ensemble 
predicted by equilibrium Statistical Mechanics. This im- 
portant issue, either for an isolated system per se or for an 
isolated system-plus-bath composite, has been recently 
addressed e.g. in 

Hi, EM 

but will not be considered 

here in any further detail. 

A first key point of our approach is the modeling con- 
cept that any given experimental system gives rise 
to a specific, well defined statistical ensemble (density 
operator) the details of which are, however, unknown 



in practice. Typically, one only knows that in the ini- 
tial state certain (usually macroscopic) observables are 
relatively sharply distributed about some approximately 
known mean values. A main challenge of the theory is to 
properly cope with this lack of knowledge. A second key 
point is not to modify or approximate in any way the ex- 
act Quantum Mechanical time evolution. The third key 
point is to focus on experimentally realistic observables, 
exhibiting a finite range and a finite resolution. 

For the rest, our present approach still covers essen- 
tially arbitrary (generic) systems, and as such is comple- 
mentary to many recent investigations of different spe- 
cific model classes, observables, and initial conditions - 
often with a main focus on the so-called eigenstate ther- 
malization hypothesis and the role of (non-)integrability 
in this context - see e.g. [3, [H, [h4T - [35| . In particular, 
our present approach is not restricted to macroscopic ob- 
servables, in contrast to e.g. [U, HE H3- The reason is 
that equilibrium Statistical Mechanics in fact also covers 
microscopic observables, e.g. the position, velocity etc. 
of one specific "tracer" particle within an fluid, and its 
pertinent (probabilistic) predictions are in perfect agree- 
ment with numerical simulations and with the rapidly 
increasing number of experiments on single molecules, 
nano-particles etc. 



II. MODEL ASSUMPTIONS 

We consider a large (macroscopic but finite), isolated 
system, modeled in terms of a Hilbert space % and a 
time-independent Hamiltonian H : T-L — > T-L of the form 



H = J2 E nPn 



(1) 



where the P n are projectors onto the eigenspaces of H 
with mutually different eigenvalues E n and multiplicities 



(h. 



Tr{P„} 



(2) 
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For the sake of simplicity, we restrict ourselves through- 
out this paper to Hilbcrt spaces % of arbitrarily large, but 
finite dimensions. Generalizations to infinite dimensions 
are readily possible along the lines of [l2j]. In particu- 
lar, the number of distinct energy eigenvalues is given by 
some finite integer Ne and the indices n in ((T|) and ([2]) 
run from 1 to Ne- 

Besides the multiplicities/degeneracies of the energy 
levels from ([2]), a related further quantity of central im- 
portance will be the maximal degeneracy of energy gaps 

g := max \{(k, I) : E k - E l = E m - E n }\ , (3) 

where 15*1 denotes the number of elements contained in 
the (finite) set 5*. In other words, g is the maximal num- 
ber of (exactly) coinciding energy differences among all 
possible pairs of distinct energy eigenvalues. In partic- 
ular, g = 1 is tantanmount to the absence of any de- 
generate energy gaps: For any given m ^ n (and thus 
E m 7^ E n ) all energy gaps Ek — Ei arc different from 
E m — E n , except for k = m and I = n. 

System states are described by density operators p(t), 
evolving according to p(t) = U t p(Q)U} with propagator 
Ut := cxp{— \Ht] and H = 1. With ([T]) we thus can 
conclude that 

= Yl <°™»(°) ex P[-i(£m ~ E n )t] , (4) 

mn 

where we have introduced the auxiliary operators 

Pmn 

(t) := P mP {t)P n . (5) 

Observables are modeled by self-adjoint operators A 
with expectation values Tr {p(t) A}. Similarly as in ([1}, 
any such A can be written in the form 

A = ^a v Q„, (6) 

V 

where Q u are the projectors onto the eigenspaces of A 
with Na mutually different eigenvalues a v . 

A single Quantum Mechanical measurement process 
of the observable A on the system in the state p(t) thus 
results in one of the Na possible outcomes a„, and the 
probability to obtain the specific outcome a„ is given by 

q v {t) := Tv{p{t)Q v } . (7) 

Upon repeating the same experiment (same measurement 
and system state) N rep times, every outcome a v is thus 
obtained about N v :— q v {t)N rep times with typical sta- 
tistical fluctuations/uncertainties of the order of y/N v . 
In other words, the "true" values of q v (t) can be deter- 
mined by means of N rep measurements up to a statistical 
uncertainty of about 



While in principle, one and the the same experiment 
can be repeated arbitrarily often, in practice (under ex- 
perimentally realistic conditions) the number of repeti- 
tions N rep must remain "reasonable" . To be on the safe 
side, we henceforth take for granted the quite weak bound 

N rep < 10 20 . (9) 

In fact, also N rep < 10 100 or N rep < 10 1000 would still 
be no problem in our later considerations, and it seems 
convincing that a theory which does not go very much 
beyond that will do. 

As pointed out e.g. in (3(| [38Tj44l | it is not necessary 
to admit any arbitrary self-adjoint operator A in order 
to model real experimental measurements. Rather, it is 
sufficient to focus on experimentally realistic observables 
in the following sense [y, [|| [45[ : Any observable A must 
represent an experimental device with a finite range of 
possible outcomes of a measurement, 

Aa := max (i/>\A\ip) - min (ip\A\ip) = a max - a mm , 
W\ =1 11-011=1 

(10) 

where the maximization and minimization is over all nor- 
malized vectors \ip) E H and where a max and a m i n are 
the largest and smallest eigenvalues of A. Moreover, this 
working range of the device must be limited to ex- 
perimentally reasonable values compared to its resolution 
limit 5A. Due to the manifold and often not very well un- 
derstood contributions to the uncertainty SA of a real ex- 
perimental measurement, quantitative estimates are not 
easy. Yet it seems reasonable to claim that uncertainties 
of less than, say 10 -20 , of the instrumental range Aa are 
impossible to achieve, i.e. 

Aa/SA < 10 20 . (11) 

Similarly as below ©, also exponents of 100 or 1000 are 
still admissible, and to go even far beyond that seems 
not necessary to faithfully model any real or numerical 
experiment. 

A particularly interesting and important example is 
a digital instrument, displaying d,A decimal places. In 
this case, we have that Na = 10 dA and we can assume 
without loss of generality that a v — SA [v — 1) with 
v = 1, Na = l0 dA and thus A A = (N A - 1) SA. It is 
quite convincing that every experimentally realistic mea- 
surement device, both digital and analoge, can be satis- 
factorily modeled along these lines. Moreover, all mea- 
surements known to the present author yield less than 20 
relevant digits, i.e. 

Aa/SA = N A - 1 < 10 20 , (12) 

in accordance with our general estimate (|11[) . 

Next we focus on the specific observable A = P n , de- 
scribing according to (fT]) the population of the (possibly 
degenerate) energy level E n with expectation value (oc- 
cupation probability) 



Sq v (t) = Jq v (t)/N, 



(8) 



Pn := Tr{P„p(t)} 



(13) 



3 



Note that P n commutes with H from ([]} and hence the 
level populations p n are conserved quantities. 

For a system with / degrees of freedom there are 
roughly 10°^ energy eigenstates with eigenvalues in ev- 
ery interval of 1J beyond the ground state energy, see e.g. 
(4a | or section 2.1 in !9']. The same estimate carries over 
to the number of energy eigenvalues under the assump- 
tion that their multiplicities ^ are much smaller than 
I0°(f). For a macroscopic system with / = 0(1O 23 ), the 
energy levels are thus unimaginably dense on any decent 
energy scale and even the most careful experimentalist 
will not be able to populate only a few of them with 
significant probabilities p n , in particular after averaging 
over the ensemble, i.e. over of many repetitions N rep of 
the experiment. In the generic case we thus can conclude 
0, [tl that even if the system's energy is fixed up to 
an extremely small experimental uncertainty, and even 
if the energy levels are populated extremely unequally, 
we still expect that even the largest ensemble-averaged 
level population p n will be extremely small (compared to 
^2 n Pn = 1) and typically satisfy the very rough estimate 

max p n = lCT° (/) . (14) 



III. DEFINITION OF EQUILIBRATION 

Generically, the statistical ensemble pit) is not station- 
ary right from the beginning, in particular for an ini- 
tial condition p(0) out of equilibrium. But if the right 
hand side of (|4]) depends on t initially, it cannot approach 
for large t any time-independent "equilibrium ensemble" 
whatsoever. In fact, any mixed state pit) returns arbi- 
trarily close (with respect to some suitable distance mea- 
sure in Hilbert space) to its initial state p(0) for certain, 
sufficiently large times t, as demonstrated for instance in 
appendix D of reference [13] . We emphasize, that these 
arbitrarily close recurrences do not refer to pure states 
only (as in the classical Poincare recurrences) but rather 
to arbitrary statistical ensembles pit). 

As an example, we focus on the simplest case that 
all eigenvalues E n in (p} are non-degenerate and thus 
every projector P n is of the form |n)(n|, where \n) is 
the (unique) eigenvector belonging to the eigenvalue E n . 
Next, we consider any p(t) which is (at least initially) not 
completely independent of t. So, according to (|4]) there 
must exists a pair of indices m, n with the property that 
(m|p(0)|n) ^ and fl := [E n - E m ]/h ^ 0. Then, by 
focusing on the specific observable 

A = B + £ f , B := \m){n\/{m\p{Q)\n) (15) 

it readily follows from @ that 

Tr{p(t)A} = 2 cos(nt) . (16) 

In other words, the ensemble p{t) exhibits permanent 
oscillations rather than equilibration, at least as far as 
the specific observable (TT51) is concerned. 



The main implication of the two previous paragraphs 
is that equilibration cannot be true (and hence cannot 
be proven) in full generality and rigor. Put differently, 
Quantum Mechanics and equilibration are strictly speak- 
ing incompatible. Equilibration can at most approxi- 
mately hold true for a restricted class of observables A 
and initial conditions p(0). We therefore focus on the 
weaker notion of equilibration from [6l-ll2l| , which merely 
requires the existence of a time-independent "equilibrium 
state" p (density operator) with the property that the 
difference 

a(t) := Ti-{p(t)A} - Ti{pA} (17) 

between the true expectation value Tr{p(t)A} and the 
equilibrium reference value Tr{pA} is unresolvably small 
for the overwhelming majority of times t contained in any 
sufficiently large (but finite) time interval [0, T). In other 
words, the expectation values Tr{p(i)A} may still exhibit 
everlasting small fluctuations around their "equilibrium 
values" Tr{/L4}, as well as very rare large excursions 
away from equilibrium (including the above mentioned 
recurrences [2|), but quantitatively these fluctuations are 
either unobservably small or exceedingly rare on any re- 
alistic time scale after initial transients have died out. 
(Note that those initial transients become irrelevant if T 
is chosen large enough.) The compatibility of the specific 
example (1151) . (|16[) with this notion of equilibration will 
be further discussed in Sect. IVII 

As we will see below, such an equilibrium ensemble p 
indeed exists under quite weak conditions, and is given 
by the following definition 

P '-=^2Pnn , (18) 

n 

where the operators p nn := p nn (t) are defined via ([5]) 
and where the time-arguments of p n n{i) have been omit- 
ted since these operators are in fact time independent 
according to (@}, ©. 

Intuitively, (|18p may be viewed as the infinite-time av- 
erage of p{t) from (U)), but also as the time- independent or 
the "diagonal" part of p(t). Yet it should be emphasized 
that the definition (fT5|) is formally self-contained and un- 
ambiguous without any such "interpretation" . In partic- 
ular, p in (|18p is well-defined, independently of whether 
the infinite-time average of pit) actually exists or not. 
Note that also physically, the "equilibrium state" p is 
mainly a theoretical "vehicle" without a direct correspon- 
dence in the real system. 

In the following, we will see that the above mentioned 
"equilibration properties" of a it) from (fl"T| indeed can be 
demonstrated by adopting this specific choice (fl"8|) of the 
"equilibrium state" p. To do so, the quantity of foremost 
interest turns out (as one might have expected) to be the 
time-averaged variance 

<a 2 (t)> T :=i^ dta 2 it) (19) 
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following from (IT7|) and ([TB]) . Considering and estimating 
such averages has a long tradition, see e.g. pi. [l4l |48M52| . 
Substantial progress along this line has been achieved 
quite recently in the works [6l-[l2j|. In the following we 
summarize and further develop this approach. 



IV. EQUILIBRATION OF MEAN VALUES 

Closely following the line of reasoning of Short and 
Farrelly [lj, we introduce Q, ([XT]), (jig)) into (H]) to 
conclude that 

(a 2 {t)) T = (| E Tr{ Pmn A}exp[-i(E m - E n )t]^ t , 

(20) 

where p mn (0) is abbreviated as p mn and the sum runs 
over the finite set of pairs of labels 

G := {(m,n) : m,n e [1,...,N E ], m f n) . (21) 

Defining for any a € Q with a — (m, n) the two quantities 

G a := E m - E n , v a := Tr{p mn A} (22) 

and introducing the self-adjoint, non-negative matrix M 
with matrix elements 



M a p := (exp [i(G a - Gp)t]) T 
we can rewrite (|2l7|) as 



(23) 



(24) 

Denoting by ||M|| the standard operator norm of the ma- 
trix M (tantamount to the largest eigenvalue of M) it 
follows that [HI 

(a 2 (t)) T < \\M\\ E K\ 2 = IIMII E \^{PmnA}\ 2 (25) 

Obviously, every matrix element (I23|) exhibits a well- 
defined limit as the averaging time T (cf. ([Ti?]) ) tends to 
infinity. Thanks to our restriction to finite dimensions be- 
low (JT|) it then follows that also the norm ||M|| converges 
for T — > oo. Hence, ||M|| can be bounded from above 
for all sufficiently large (but finite) T. Quantitatively, 
for non-degenerate energy gaps, i.e. q = 1 in (|3|), one 
can infer from (|23p that M approaches the unit matrix 
for T — > oo and hence ||M|| 1 = g. As demonstrated 
in in [IJ [l2j], the same result ||M|| — > g for T — > oo re- 
mains valid for arbitrary g. We thus can conclude that 
||M|| < 2g for all sufficiently large T, where g is the 
maximal degeneracy of energy levels from ([3]). Observ- 
ing that the sum in (|25|) can only grow by extending it 
to all m, n € [1, . . . , Ne] it follows that 



(a 2 (t)) T <2g]T|Tr{p„ m A}| 2 



(26) 



for all sufficiently large T. We remark that this is the only 
place where our restriction to finite dimensions below ([T]) 
is essential. While it is quite plausible that everything 
still "goes well" in the limit of infinite dimensions [f| Q , 
a more detailed justification is actually not so obvious 
and has been worked out in 12]. 

Next we employ (J5J) for p mn — p mn (0), the property 
= P m of any projector P m , and the cyclic invariance 
of the trace to conclude that 



Tr{p mn A} = Tr{P m p(0)P n AP m } 



(27) 



Since p := p(0) is a self-adjoint, non-negative operator, 
one can infer that there exists a self-adjoint, non-negative 
operator, which we denote by p 1 ^ 2 , and which satisfies the 
relation p 1 / 2 p 1 / 2 = p. Observing that 

P mP P n AP m = (P ro/ 9 1 / 2 )(p 1 /2p nj4 p ro ) (28 ) 

and exploiting the Cauchy-Schwarz inequality 

\Tr{B^C}\ 2 < Tr{B i <B}Tr{C i <C} (29) 

for the scalar product Tr{B^C} of arbitrary operators B 
and C it follows with ([5]), ([13]), and <j27j) that 

|Tr{ Pmn A}| 2 < Pm Tr{p nn AP m A} (30) 

and hence 

\Tv{p mn A}\ 2 < Tr{(J2pnn)A(£ > Pn l P m )A}(Zl) 



E v <* exp [-iG a t] ) = E v^MapVfi . Introducing 



(32) 



and taking into account (|18[) we thus obtain 

E |Tr{p m „A}| 2 < Tv{pAujA} = Tr{ujApA} . (33) 



Finally, we observe that 

Tr{BC} < ||S||Tr{C} 



(34) 



for arbitrary self-adjoint, non-negative operators B and 
C, where 



max{(^|B|^) : {^W) = 1} 



(35) 



denotes the standard operator norm (largest eigenvalue 
of B). Since both ui and ApA are self-adjoint and non- 
negative, it follows that 



Tr{uApA} < \\u\\ Tr{ ApA} = \\u\\ Tr{pA 2 } 



(36) 



The first term on the right hand side equals max n p„ 
according to (|32[) . In view of (f26| and (|33|) we thus can 
conclude that 



(<J 2 {t)) T < 2Tr{pA 2 }g maxp„ 



(37) 
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for all sufficiently large T. Since both p and A 2 are self- 
adjoint and non-negative we can exploit (|34[) to conclude 
that Tr{p,4 2 } < \\A 2 \\ Tr{p}. With Tr{p} = 1 it follows 
that Ti{pA 2 } < \\A 2 \\. Next we note that a(t) in ((T7J) and 
thus (cr 2 (i)) T in (fT9|) remain unchanged if A is replaced 
by A + c 1 , where c is an arbitrary real number and 1 
the identity operator. Choosing c so that a m i n — —a max 
and thus a max = Aa/2 in (|10[) it follows with (1351) that 
j|^4j| = Aa/2. Altogether, we finally obtain the result 
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{° 2 {t)) T < - A\g maxp„ 



(38) 



N rep times, the "true" frequency q^N rep of the outcome 
a„ can thus be determined up to a statistical uncertainty 
of about 



5q v ■= \ q v /N, 



v rep 



Focusing on the special choice e 
from (1351)- (Hill that 



(44) 

8q u we can conclude 



T 



<2N r 



, g max p T , 



(45) 



for all sufficiently large T. 

Next we define for any given 5A > and T > the 
quantity 

T SA := |{0 < i < T : \Tr{p{t)A} - Tv{pA}\ > 5A}\ , 

(39) 

where |5| denotes the size (Lebesgue measure) of the set 
S. According to (1T7|) . Tsa thus denotes the measure of 
all times t £ [0,T] for which \a(t)\ > SA holds true. It 
follows that a 2 (t) > SA 2 for a set of times t of measure Tsa 
and <J 2 (t) > for all remaining times t in [0,T]. Hence 
the temporal averag e flTjJJ) must be at least 5A 2 T SA /T and 
we can conclude from (1581 that 



Tsa < 1 
T ~ 2 



&4 



5 maxp„ 



(40) 



for all sufficiently large T. This is the first main result 
of our present paper. A more detailed discussion of its 
physical content is postponed to Section |VT1 



for all sufficiently large T. According to fill) , the left 
hand side of (j45|) quantifies the fraction of all times t € 



[0, T] for which one may find by means of N rep measure- 
ments a statistically significant difference between the 
two frequencies q v {t) — Trjp^Qj,} and q v = Tr{pQ u } of 
observing the outcome a v . For all other times t € [0, T], 
the "true" frequency q v {t) cannot be distinguished from 
the one theoretically predicted by means of p. Denoting 
by T* the measure of all times t £ [0, T] for which such a 
distinction is possible at least for one of the Na different 
indices v, one readily sees that T* < T$q v>v and thus 
with Ell) that 



T 



< 2 N rep Na g maxp rt 



(46) 



for all sufficiently large T. For all other time t £ [0, T), it 
is impossible to observe by means N rep measurements of 
A any kind of statistically significant difference between 
p(t) and p. 



V. EQUILIBRATION OF INDIVIDUAL 
MEASUREMENT OUTCOMES 

Next we focus on the special case that A happens to 
coincide with one of the projectors Q v in Similarly 
as in (|39|) . we define for any given e > and T > the 
quantity 

r e ,„ := |{0 < t < T : \Tt{p(t)Q v ] - Tr{pQA\ > c}| , 

(41) 

and similarly as in (|40p we can conclude from (|37[) and 
Ql = Qv that 

T 2 

~Tf~ < -o Tr{pQ„} g maxp„ (42) 

for all sufficiently large T. 

Returning to general observables A of the form ^ we 
denote, similarly as in (J7]), by 

q v := Tt{pQu} (43) 

the probability that a single measurement of A would 
result in the outcome a v if the system were in the state p. 
Similarly as in ([8J, upon repeating the same experiment 



VI. DISCUSSION 

The first main result of our paper is the relation (1401) . 
In fact, the same result has been obtained already be- 
fore by Short and Farrelly in [111, however by means of a 
quite different line of reasoning [54j. In particular, these 
authors focus in a first step solely on pure states and only 
in the end extend their result to mixed states via purifi- 
cation. On the other hand, our pre sent line of reasoning 
is very similar to the one from [12|], but the final result 
is significantly stronger. 

According to (|3"9"|) . the left hand side of (j40l) repre- 
sents the fraction of all times t £ [0, T) for which there is 
an experimentally resolvable difference between the true 
expectation value Tr{p(t)A} and the time- independent 
equilibrium expectation value Tr-fjoA}. On the right hand 
side of (|40|) . Aa/SA is the range-to-resolution ratio, which 
may be very large but still remains bounded according 
to (JXT]) for experimentally realistic observables. The next 
factor, g, is the maximal degeneracy of energy gaps from 
([3]). Finally, max„p n represents the largest, ensemble- 
averaged occupation probability of the (possibly degen- 
erate) energy eigenvalues E n , see (JT3J) . Typically, one 
expects that the rough upper bound (j!4[) applies, except 
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if certain energy eigenvalues are so extremely highly de- 
generate that the multiplicities defined in ^ severely re- 
duce the pertinent energy level density compared to the 
non-degenerate case, see above (fH|) . 

For a system with sufficiently many degrees of free- 
dom / and no exceedingly large degeneracies of energy 
eigenvalues and energy gaps 55|, we thus can conclude 
from (|4T))) with (fTI]) and (fT4")l that the system behaves 
with respect to any experimentally measurable expecta- 
tion value exactly as if it were in the equilibrium state p 
for the overwhelming majority of times within any suffi- 
ciently large (but finite) time interval [0,T]. In particu- 
lar, T must obviously be much larger than the relaxation 
time in case of a far-from-equilibrium initial condition 
P (0). 

Next we turn to the second main result of our pa- 
per, namely Eq. (146p . Its left hand side represents 
the fraction of all times t £ [0, T] for which one may 
find by means of N rep measurements of A a statisti- 
cally significant difference between the "true" probabil- 
ities q v (t) = Tr{p(t)Q v } and the corresponding "equi- 
librium probabilities" q v = Tr{/?Q„} of observing any of 
the Na different possible outcomes a v of the Quantum 
Mechanical measurement process. According to ([9]) , (TT21) . 
and (fnj we can conclude from (j4U)) that for a system with 
sufficiently many degrees of freedom / and no exceed- 
ingly large degeneracies of energy eigenvalues and energy 
gaps, the true state p(t) is practically indistinguishable 
from the equilibrium state p after initial transients have 
died out. More precisely, since a real experiment can 
only be repeated a "reasonable" number of times (see 
©), the observable differences between p(t) and p are 
either smaller than the remnant statistical uncertainties 
or so unimaginably rare in time that they never actually 
occur in practice. 

The main difference between (|40|) and (|46|) is that the 
former refers to the expectation value of A, i.e. the mean 
value over all the actually observed outcomes a v of the 
single repetitions of the experiment. In particular, statis- 
tical uncertainties due to finite numbers of experiments 
are not taken into account, i.e. N rep — > oo is implicitly 
taken for granted. The result (j4"fj)) goes substantially be- 
yond (|40|) in two respects. First, it properly accounts for 
a finite number of experiments. Second and even more 
importantly, not only the first moment (mean value) but 
actually the full probability distribution of the different 
possible measurement outcomes is covered by the result 

(USD. 

It is instructive to illustrate the above general consid- 
erations by means of our specific example from (|15p , (| 1 C5[) . 



Recalling that in this example we are dealing with a non- 
degenerate Hamiltonian |T]) with eigenvectors \n), one 
readily sees that the spectrum of A from (fT5|) consists of 
the two (non-degenerate) eigenvalues a± = ±\p nm (0)\~ 1 , 
and one (highly degenerate) eigenvalue eto = 0. In other 
words, we have Na = 3 distinct eigenvalues in ©. Fur- 
thermore, one can infer, similarly as in (1281) and exploit- 
ing d2U), that |(m|p(0)|n)| 2 < (m\p(0)\m) (n\p{0)\n). Ac- 
cording to (fT"3|) . the latter product equals p m p n - It fol- 
lows that \a±\ > l/max n p„ and hence with (|f Op that 
Aa > 2/max n p„. For experimentally realistic initial 
conditions we can infer with (fT4")) that > 0(10?). On 
the other hand, in order to resolve the oscillations from 
([in]), the experimental uncertainty 5A must not exceed 
unity. 

In other words, while it is clearly impossible to deny 
the existence of the permanently oscillating expecta- 
tion values (HHJ), they cannot be actually "seen" without 
admitting unrealistically large range-to-resolution ratios 
Aa/SA or statistical ensembles with unrealistically large 
energy level populations max n p n . 

The same thing may alternatively also be viewed as 
follows: Any single (ideal) measurement process always 
results in one of the three outcomes a+, a_, or ao- Recall- 
ing that a + = — a_ > l/max„p„ and ao = it follows 
that an infeasible number of measurements N rep would 
be needed to resolve the order-one variations of the en- 
semble average (fTtJ)) . 

The above example also suggests that our assumptions 
of experimentally realistic observables, number of mea- 
surements, and initial conditions (or some similar restric- 
tions) are unavoidable for taming the oscillations in (0| 
and thus overcoming the fundamental incompatibility of 
basic Quantum Mechanics with rigorous equilibration, as 
detailed in Sect. IIII1 In other words, it seems not an ar- 
tifact of our present (or of any other) specific line of rea- 
soning but rather a feature of "physical reality" (within 
the realm of standard Quantum Mechanics) that invok- 
ing "our human inability" is unavoidable to demonstrate 
and understand equilibration. 



Acknowledgment 

Collaboration with Michael Kastner on closely related 
issues is gratefully acknowledged. This work was sup- 
ported by Deutsche Forschungsgemeinschaft under grant 
RE1344/7-1. 



[1] Sklar L 1993 Physics and Chance, Cambridge University 
Press 

[2] Penrose O 1979 Rep. Prog. Phys. 42 129 
[3] Farquhar I E 1961 Ergodic Theory in Statistical Me- 
chanics, Interscience, NY; Jancel R 1969 Foundations of 



Classical and Quantum Statistical Mechanics, Pergamon, 
London 

[4] Focus Issue Dynamics and thermalization in isolated 
quantum many-body systems, edited by Cazalilla M A 
and Rigol M 2010 New J. Phys. 12 



7 



[5] Polkovnikov A, Sengupta K, Silva A and Vengalattore M [32 
2011 Rev. Mod. Phys. 83 863 [33 

[6] Reimann P 2008 Phys. Rev. Lett. 101 190403 

[7] Linden N, Popescu S, Short A J and Winter A 2009 Phys. [34 
Rev. E 79 061103 

[8] Linden N, Popescu S, Short A J and Winter A 2010 New [35 
J. Phys. 12 055021 

[9] Reimann P 2010 New J. Phys. 12 055027 [36 
[10] Short A J 2011 New J. Phys. 13 053009 

[11] Short A J and Farrelly T C 2012 New J. Phys. 14 013063 [37 
[12] Reimann P and Kastner M 2012 New J. Phys. 14 043020 
[13] Goldstein S, Lebowitz J L, Mastrodonato C, Tumulka 

R and Zanghi N 2010 Phys. Rev. E 81 011109; Gold- [38 
stein S, Lebowitz J L, Tumulka R and Zanghi 2010 Eur. [39 
Phys. J. H 35 173; Goldstein S, Lebowitz J L, Mas- 
trodonato C, Tumulka R and Zanghi N 2010 Proc. R. Soc. [40 
London Ser. A 466 3203; Tasaki H 2010 The approach [41 
to thermal equilibrium and "thermodynamic normality" 
Preprint laFXiv:1003.5424V 4 [42 
[14] Tasaki H 1998 Phys. Rev. Lett. 80 1373 
[15] Srednicki M 1999 J. Phys. A : Math. Gen. 32 1163 
[16] Cazalilla M A 2006 Phys. Rev. Lett. 97 156403 [43 
[17] Rigol M, Dunjko V, Yurovsky V and Olshanii M 2007 [44 

Phys. Rev. Lett. 98 050405 
[18] Kollath C, Lauchli A M and Altman E 2007 Phys. Rev. [45 
Lett. 98 180601 [46 
[19] Manmana S R, Wessel S, Noack R M and Muramatsu A 

2007 Phys. Rev. Lett. 98 210405 
[20] Cramer M, Dawson C M, Eisert J and Osborne T J 2008 [47 

Phys. Rev. Lett. 100 030602 
[21] Barthel T and Schollwock U 2008 Phys. Rev. Lett. 100 [ 

100601 [49 
[22] Rigol M, Dunjko V and Olshanii M 2008 Nature 452 854 [50 
[23] Roux G 2010 Phys. Rev. A 81, 053604 [51 
[24] Pal A and Huse D A 2010 Phys. Rev. B 82 174411 [52 
[25] Biroli G, Kollath C and Lauchli A M 2010 Phys. Rev. 

Lett. 105, 250401 [53 
[26] Gogolin C, Miiller M P and Eisert J 2011 Phys. Rev. 

Lett. 106 040401 [54 
[27] Ikeda T N, Watanabe Y and Ueda M 2011 Phys. Rev. E 

84 021130 [55 
[28] Bafiuls M C, Cirac J I and Hastings M B 2011 Phys. Rev. 

Lett. 106 050405 
[29] Campos Venuti L, Jacobson N T, Santra S and Zanardi 

P 2011 Phys. Rev. Lett. 107 010403 
[30] Jacobson N T, Campos Venuti L and Zanardi P 2011 

Phys. Rev. A 84 022115 
[31] Kastner M 2011 Phys. Rev. Lett. 106 130601 



Ji K and Fine B V 2011 Phys. Rev. Lett. 107, 050401 
Canovi E, Rossini D, Fazio R, Santoro G E and Silva A 
2011 Phys. Rev. B 83, 094431 

Kastner M 2012 Cent. Eur. J. Phys. 

DOI:10.2478/sll534-011-0122-4 (in press) 

Rigol M and Srednicki M 2012 Phys. Rev. Lett. 108, 

110601 

Khinchin A Y 1960 Mathematical Foundations of Quan- 
tum Statistics (New York: Graylock Press). 
Wightman A S 1971, pp 1-32 in Statistical Mechanics at 
the turn of the Decade, edited by E. G. D. Cohen, Marcel 
Dekker Inc., New York 
Wigner E P 1963 Am. J. Phys. 31 6 
van Kampen N G 1992 Stochastic Processes in Physics 
and Chemistry (Amsterdam: Elsevier) chapter XVII. 7 
Sugita A 2007 Nonlinear Phenom. Complex Syst. 10 192 
Penrose O 1970 Foundations of Statistical Mechanics 
(Oxford: Pergamon) chapter 1 

Martin-L6f A 1979 Statistical Mechanics and the Founda- 
tions of Thermodynamics (Lecture Notes in Physics vol 
101) (Berlin: Springer) 
Georgii H-O 1995 J. Stat. Phys. 80 1341 
Popescu S, Short A J and Winter A 2006 Nature Phys. 
2 754 

Reimann P 2008 J. Stat. Phys. 132 921 
Landau and Lifshitz, Statistical Physics, Pergamon, Ox- 
ford 1980; Diu B, Guthmann C, Lederer D and Roulet B 
1996 Elements de Physique Statistique (Paris: Hermann). 
Hobson A 1971 Concepts in Statistical Mechanics, Gor- 
don and Breach, New York 

Ludwig G 1958 Z. Phys. 150, 346 (1958); 152, 98 
Bocchieri P and Loinger A 1959 Phys. Rev. 114, 948 
Peres A 1984 Phys. Rev. A 30 504 
Deutsch J M 1991 Phys. Rev. A 43, 2046 
Srednicki M 1996 J. Phys. A: Math. Gen 29, L75; Sred- 
nicki M 1999 J. Phys. A: Math. Gen 32, 1163 
Goldstein S, Lebowitz J L, Tumulka R and Zanghi N 
2006 J. Stat. Phys. 125 1197 

In fact, the estimate obtained in [lj| is even somewhat 
tighter than that in (|40|) . 

Note that no such degeneracies are actually encountered 
in the generic case @, 0, 0, Q3, 0, , i.e. in the ab- 

sence of special reasons like additional conserved quanti- 
ties (e.g. due to (perfect) symmetries or when the system 
consists of non-interacting subsystems) or fine-tuning of 
parameters ("accidental degeneracies"). 



